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Abstract 

The diffeomorphism symmetry of general relativity leads in the canonical formulation to con- 
straints, which encode the dynamics of the theory. These constraints satisfy a complicated algebra, 
known as Dirac's hypersurface deformation algebra. This algebra has been a long standing challenge 
for quantization. One reason is that discretizations, on which many quantum gravity approaches rely, 
generically break diffeomorphism symmetry. 

In this work we find a representation for the Dirac constraint algebra of hypersurface deformations 
in a formulation of discrete 3D gravity and for the flat as well as homogeneously curved sector of 
^vj ' discrete 4D gravity. In these cases diffeomorphism symmetry can be preserved. Furthermore we 

£*\j , present different versions of the hypersurface deformation algebra for the boundary of a simplex in 

arbitrary dimensions. 
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i 1 Introduction 

5-1 ' 

Diffeomorphism symmetry is the fundamental gauge symmetry of general relativity, deeply entangled 
with its dynamics. In the canonical framework this gauge symmetry leads to the Hamiltonian and diffeo- 
morphism constraints. They generate multi-fingered time evolution - that is the equal time hypersurface 
can be deformed and pushed forward in various ways. These deformations of the hypersurface satisfy 
a geometric algebra which is reflected in the Poisson algebra of the constraints [I], known as Dirac's 
Q\ ' hypersurface deformation algebra [2]. 

l/~) . An important goal for the canonical quantization program is to find a quantum representation of 

,—4- ' the algebra [31 H] [5]. It has not been reached yet, even in the case of 3D gravity which is a topological 

theory and is well understood in many of its aspects. The Dirac algebra is universal [T], i.e. it is the 

same algebra for any theory of hypersurfaces embedded in a higher dimensional manifold. To understand 

possible quantum representations of this algebra would therefore be beneficial not only for gravity. 

One reason that such a representation is still missing, is that many approaches use discretizations 
either from the outset or as auxiliary structures to e.g. regularize the quantum constraints. However 
discretizations generically break diffeomorphism symmetry [5J |S] , even classically and in simple systems 
[TTJl [TTj . For instance, the canonical analysis of Regge gravity [H] shows that in this case constraints are 
replaced by proper evolution equations [13] . In 3D, Regge gravity [14] does however provide a symmetry 
preserving discretization, which has to be adjusted to the presence of the cosmological constant [15] 
though. Thus we might hope that a quantum representation of the constraint algebra might be obtainable 
at least in this case. Indeed quantizations of the constraints are available [THl [TTJ [T51 [TH] , although a 
representation of the Dirac algebra, even classically, has not been found so far. 

3D quantum gravity [201 (HI HH H3] had been used as a simplified model for 4D gravity for some time. 
It can be formulated in a number of ways, in traditional (canonical) ADM form, for instance [23] , in 
Chern-Simons form [551 HI] or as a so-called BF theory, which coincides with the Palatini formulation 
of 3D gravity. 

We will work with the latter version, as it is nearest to the variables and methods used in loop 
quantum gravity, see for instance [25[ . There has been also a considerable amount of work on the loop 
quantum gravity version of 3D gravity [26l [TT1 [2TJ [28j [19] . Most of this work, with the exception of [17] 
makes use of the fact that the Hamiltonian and diffeomorphism constraints can be recast into the flatness 
constraints. Whereas the first set satisfies the Dirac algebra of hypersurface deformations, the second set 
is Abelian. 

Despite a good understanding of discrete 3D gravity there is so far no complete discussion or defini- 
tion of Dirac's hypersurface deformation algebra for discrete surfaces available. The works [3H1 [30] are 



incomplete for various reasons, with [29 arguing (incorrectly) that such an algebra has to be necessarily 
non-local. 

Thus we will provide here an explicit classical realization of the hypersurface deformation algebra for 
discrete 3D gravity. We discuss the geometric interpretation of the constraints and use these insights to 
define (various) hypersurface deformation algebras also for discrete 4D gravity (and higher dimensional 
discretized surfaces). In this case we have however to restrict to the flat sector (or homogeneously curved 
sector if a cosmological constant is present) of these theories, as otherwise diffeomorphism symmetry will 
be broken by the discretization. Nevertheless this justifies the hope that some results of the 3D theory 
can be applied to the 4D case. This might even hold for the quantum theory [31] . 

The organization of the article is as follows. In the Section [2 we review the continuum formulation 
of 3D gravity, in particular the Dirac algebra, and we present in the Section [3] the discretization we are 
going to use, inherited from loop quantum gravity. In the Section 0] we introduce the Hamitlonian and 
diffeomorphism constraints, using a natural space-time splits coming from the loop quantum gravity 
fluxes. We discuss the form of the algebra and the main tool for its computation in the Section \5\ In 
the Sections [6] and we give the explicit Poisson brackets between the constraints associated to a single 
3-valent vertex and to a face respectively. The hypersurface deformation algebra is written in terms of 
lengths and dihedral angles in the Section HI The Section [9] shows that the Thiemann trick, used to dchnc 
the quantum Hamiltonian in loop quantum gravity (in particular (17) is also available in the discrete 
setting, and can be used to write our Hamiltonian constraints. We then move on to higher dimensions, 
starting with the 4D case in the Section [TOl in which we argue that a closed algebra is available in the 
flat sector of the theory. We introduce the corresponding constraints in the Section El which work in 
arbitrary dimension on the boundary of a d-simplex. The constraint algebra is found in the Section 1121 
and generalized to the case of a non-vanishing cosmological constant (describing homogeneously curved 
space) in the Section El We close the paper with a discussion in the Section El The Appendix [2] 
provides details on the evaluation of the Poisson brackets for 3D gravity, relevant for the Sections [6] and 
[3 and the Appendix [B] contains background material on the afhne metric in a simplex that we use in the 
Sections El El 

2 Continuum formulation of 3D gravity 

3D (Euclidean) gravity can be formulated in first order form, in which the co-triad efj, k — 1,2,3 and a 
connection one-form A k are the basic variables. Here \x, p, . . . denote space time indices and i, k, j internal 
(su(2) algebra) indices, which are raised and lowered with the internal Euclidean metric Ski- With the 
curvature tensor 

F^ = d,Ai~d u Al + e l ^A M A uk (2.1) 

the action for a 3D manifold M = SxM can be written as 

S=\ [ e al Fl v e^d 3 x = f Ud Q A aj € ah + e 0j ^ b e ab + A j Q {d a e hj + e j i m A l a e b n )e ab ) d^dx . 

(2.2) 
Here e ab , e CT/iI/ are totally antisymmetric Levi-Civita tensor densities and e J is the Levi-Civita tensor in 
the internal indices. We use the Einstein summation convention. In the second line in (J2.2I) we isolated 
the time components // = from the spatial components (i — a,b, describing the spatial hypersurface £ 
by a: = const.. 

Thus we find the momenta 

E?(x t t) := ■ = i ab e bl (x,t) (2.3) 

giving the canonically conjugated pair 

{Ai(x),E b k (y)} = S{5 b a 6(x,y) . (2.4) 

The other Poisson brackets vanish: {A k (x),A J b (y)} = and {E^(x), E b (y)} = 0. 



The time components of e, A are Lagrangian multipliers for the action (|2.2p . varying this action with 
respect to e J and A^ we find the GauB and flatness constraints 

Qi := \T ab] e ab = d a E] + e jlm A l a E am =0 

& ■= KV & = £ ab (d a Al + ^ kl A k a A b ) I (2.5) 

where T l — d^ei — d v e' l u + €■> Apjevk ~ e A^e^k is the torsion of the connection A and triad e. 
These constraints are first class forming the (Galilean symmetry) algebra 

{G[A'],g[A]} = G[[A',A}} 

{T[N],gm = ?m,A}} 

{T[N'],F[N]} = (2.6) 

where G[A] = JA^Gj^x and F[N] = fN j Tj(Px. The bracket [A,B] k = e iik A l Bi is the Lie algebra 
bracket of su(2). 

The flatness constraints generate translations in the triad variables E whereas the Gaufi constraints 
generate rotations in the internal space. Gauge symmetries corresponding to diffeomorphisms arise as 
combinations of the translations and rotations, see for instance [17l [32l [25] . 

This leads to the following constraint generating spatial diffeomorphisms 



H a = -ei^ - AiG, = E]Fl b - AiGj (2.7) 



and the Hamiltonian constraint 



H = -n j ^ = -±±= e i kl E%EfF abj . (2.. 



Here r\? is the normal to the hypersurface x° — const. 

1 1 



n 



l ah ^ l E a „E\ . (2.9) 
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This form of the constraints is the same as the form of the constraints in 4D expressed in (self-dual) 
Ashtekar variables [33] . 

The term in (|2.7[) proportional to the Gaufi constraints is often omitted. Again we introduce the 
smeared constraints 

H[N}= J NHd 2 x, H[N}= I N a H a d 2 x . (2.10) 

is is 

The following Poisson bracket algebra holds modulo terms proportional to Gaufi constraints (see for 
instance [331 Hi] for the computation) 

{H[N],H[M]\ = -H[C a N] 
{H[N],H[M}} = -H[C$N\ 
{H[N],H[M}} = +H[V] with V a = q ab (Md b N - Nd b M) . (2.11) 

Here q ab is the inverse of the two-metric q ab — e^ebk, and C is the Lie derivative. 

This Poisson algebra is known as Dirac's hypersurface deformation algebra and describes the commu- 
tator of normal and tangential deformations of a hypersurface. It is universal, i.e. holds in all dimensions 
and with arbitrary field content. A special feature of the Dirac algebra is the appearance of the structure 
function q ab in the brackets between two Hamiltonian constraints. 

3 Discretization of 3D gravity 

Here we will review shortly the standard choice of discretization for 3D gravity in first order form, 

[261 EH El]. 



The internal index k in A k , e k transforms according to the fundamental representation of SO (3), i.e. 
the spin one representation. This agrees with the adjoint representation on the Lie algebra su(2), which 
suggest to introduce the Lie algebra valued forms A a = A k T k and e& = e^T k , where T k is a basis of 
su(2). We will work with T^ = T k = —jOk where Ok are the Pauli matrices!^ 

The parallel transport along a curve 7 of a Lie algebra valued object V is defined by V(s) — 
h(s)V(0)h(s)~ 1 where h(s) is the holonomy of the connection A given by 

h 7 (s) = Pexp - / A 

t/'V 

■■= V(-l)" / ds„ / dan-x---/ ds 1 A(s n )A(s n -i)---A(s 2 )A(s 1 ) (3.1) 

„_n JO JO JO 
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with A(s) = r y a (s)A a ('j(s)). A dot denotes differentiation with respect to the curve parameter s. Note 
that with this definition we have for the composition of edges e\ o e 2 the relation h ei0e2 — h e2 h ei . 

Our choice of discretization will replace the connection by holonomies g e := h e (l) along oriented edges 
e. The triad variables will be encoded into the fluxes 

E e = f h e>e *(s)e a (e*(a))(e*) a (s)(h e , e *(a)r 1 ds . (3.2) 



where e* is an edge cutting e transversally and such that (e, e*) are positively oriented. Here h e ^ e *(s) 
parallel transports vectors from the point s on e* to the source vertex v s = e(0) of e. We define h e , e *(s) 
as follows: Choose the parametrization of the curves e and e* such that the intersection point in both 
cases corresponds to s = |. Then 

h e ,e-(s) = (h e (^))- 1 -h e ,(l, S ) (3.3) 

where h e *(h,s) is the parallel transport from the parameter s to the parameter ^ along e*. 

We inserted the holonomies to ensure a local transformation behaviour of the fluxes under internal 
rotations E e — > g(v s ) E 1 ^(ii s ) _1 where g(-) denotes the gauge field on the spatial hypersurface E. 

Under the inversion of the orientation of the edge e — > e :~ e™ 1 the holonomies and fluxes transform 

as 

h s = (he)- 1 , E s = -h.E^)- 1 . (3.4) 

The Poisson brackets between holonomies and fluxes can be computed and result [341 125) in 



{(g e )MN, (g e ')M'N'} — , 

{E J e ,g e >} = Se^geT^ 

{ElE k } = bee«? k iE\ . (3.5) 

We can deduce from (13.5[) for the bracket with an inverse group element {E( ,g~ x } — —T^g' 1 as well 
as the relation {E k , E l e } = 0. 

The Gaufi constraints are assigned to vertices and are represented as 

g v = Yl E * + J2 E * ' ( 3 - 6 ) 

e:v a (e)—v e:vt(e)—v 

i.e. as sum over the fluxes associated to the outgoing (e : v s (e) — v) and the incoming e : «t(e) — v edges 
at v. The Gaufi constraints generate gauge transformations on the internal index, i.e. internal rotations. 
Later-on, in the Section[Sl we will discuss the reduced phase space with respect to the Gaufi constraints. 
From the construction of the holonomy and flux variables we can infer the following interpretation. 
The flux E e gives the components of the edge e* in a frame associated to the source vertex v s (e) of e. 
The holonomy variables g e are used to transport between the frames associated to the source vertex v s (e) 
and target vertex Vt(e) of e. Thus the Gaufi constraints impose that the fluxes around a vertex sum to 
zero, which geometrically imposes that the dual edges e* form a closed curve. For a three-valent vertex 



1 These satisfy o^o^ = Sj^l + itj^^u from which T{Tj = — i(5jjl+ o e ij k Tk an d TiTjT m = -r (<5;„ 



£<Eij m l follows. 



this will define a closed triangle. (For higher-valent vertices the edges e* making up the piecewise linear 
curve close, but will in general not span a 2-dimensional subspace.) Due to this interpretation we will 
assume that the edges and vertices are the elements of a 2-complex, which discretizes the 2D surface E. 
This latter fact means that we can also identify faces /. The 2-complex is usually assumed to arise as 
the dual of a triangulation (in which case the vertices are three-valent) or more generally a discretization 
with polygonal cells (which allows higher valent vertices). 

The flatness constraints impose that the parallel transport around (contractible) loops is trivial. The 
most elementary loops are given by the boundaries of the faces. We have furthermore to specify the 
frame, i.e. vertex, in which the holonomy is expressed. Thus we define the flatness constraints^ as 

Ff,v := \{h~ f l ~ h fv ) (3.7) 

with hf v the holonomy around a face (dual to a vertex in the triangulation) starting with the edge e C / 
which we assume to have v as a source vertex v — v s (e), i.e. hf v — • • • 9e"9e'ge- In the following, we will 
often parametrize SU(2) matrices as h — hi + h l T l , where hi — —2tr(hT l ) is the projection of h onto 
the generator T\ The inverse of h writes h~ l = h°I — h % T % . Therefore, 

F f . v = -h) v T\ (3.8) 

As can be easily checked the algebra of flatness and Gaufi constraints is again first class, with the 
flatness constraints forming an Abelian subalgebra. In the following subsection we will find combinations 
of these constraints, that can serve as discrete versions of the Hamiltonian and diffeomorphism constraints. 

4 Geometry and constraints associated to a face vertex pair 

The action of the flatness constraints on the flux associated to e is given by 

{E*,F ffV } = -\(T k hf v + h fv T k ) ~ -T k (4.1) 

with the last equation holding on the constraint hypersurface. Again we assume that hf, v = ■ ■ ■ g e "ge'ge- 
The action (|4.1j) of the flatness constraints generalizes to the fluxes E e i with e' in the boundary of /, 
if these are transported to the reference frame of v s (e). Thus Ff, v generates translations in the (to v 
parallel transported) components E k of the fluxes associated to the edges e' in the boundary of the face 
/. This face / is dual to a vertex /* in the dual triangulation, and we assume that the face is such that 
e*, the dual edge to e, points towards the dual vertex /*. 
Now contract Ff tV with a vector — n k , 

- F f , v -n:=-J2 F f,y ■■= - E(- 2 ) M-F/.^n') . (4.2) 

l k 

This defines a combination of constraints which according to 

{E k , -F f . v ■ n} ~ 2(-2) tr(T k T l n l ) = n k (4.3) 

I 

generates the translation of E k by the vector n k . 

Thus Hamiltonian (and diffeomorphism) constraints may be defined by contracting the holonomy 
around each face with some choice of normals (and vectors tangent to the hypersurface). However, 
the holonomy hf v probes the curvature around a vertex of the triangulation and there is no obvious 
notion of tangent and normal vectors to the vertex. One could proceed with some averaging over the 
normals of the adjacent triangles. That would require to specify the details of the averaging. The loop 
quantum gravity version [17] proceeds in another way and associates one Hamiltonian constraint (the 
diffeomorphisms being treated with a different method) to each dual vertex (corresponding to a triangle in 
the triangulation) . Thus there is an averaging over the three vertices of the triangle and correspondingly 
of the three holonomies involved. Note that the counting of constraints might not necessarily match up in 



2 We choose a symmetrized form of the constraints, as this is the form that is usually employed in the Hamiltonian 
constraints. It however also allows for the solution hf v = —I. One can also use Ff v = I — hf v - this will not change the 
algebra of the constraints. In the following we will understand as constraint hypersurface the subspace defined by hf v = I. 



rif v = Ei x E 2 




Figure 1: The vertex v is dual to the triangle (in solid lines) and the face / (dashed lines) is dual to a 
vertex of the triangle. Ei,E 2 are associated to the two edges meeting at v in the boudary of /. They 
describe the embedding of the triangle in R 3 , with normal iif v — {E\ x E<z)/\E\ x E 2 \. 



this case, that is one might obtain a redundant set or an incomplete set. For 4D gravity this discrepancy 
has been noted in |35j . here we see that one should be aware of possible redundancies between the 
constraints. 

Instead of averaging, we will work with an over-complete basis of constraints. Possible averagings 
can be performed afterwards and the constraint algebra adjusted, which in particular is straightforward 
if the averaging coefficients are constants. 

A vertex of the triangulation is shared by at least three triangles. Each of them has a local frame, a 
basis of which being provided by the two edge vectors of the triangle meeting at the vertex, say Ei, E 2 , 
which span the tangential directions to the triangle, and an obvious choice of normal, 



n fv = {E 1 x E 2 )/\Ei x E 2 \. 



(4.4) 



In the dual graph, the edge vectors E\,E 2 are equivalently identified by a pair face- vertex (f,v). The 
face / is dual to the vertex of the triangulation, and the vertex v labels the triangle that provides a local 
basis. Ei,E 2 are the vectors associated to the two dual edges meeting at v in the boundary of /. In the 
definition of the normal (|4.4p . we have assumed that Ei,E 2 are both outgoing (or both ingoing) at v. 
The situation is summarized in the Figure [TJ 

We can now define a Hamiltonian and diffeomorphism constraints by contracting the holonomy, or 
rather the flatness constraint Ff, v with the appropriate vectors for each pair (/, v), 



H fv 
Dl v 



-Ff v ■ nj v 



-Ff v ■ E Q 



n k h) v 



E k a h% 



{Ei x E 2 ) 



\Ex x E 2 
E k (-2)tr(h fv T k ) 



(-2)tr(h fv T k ), 



(4.5) 



Here and in the following we will sum over pairs of internal indices, even if these are not in the usual 

summation convention position (as the internal indices are raised and lowered with the trivial metric Ski)- 

To describe the algebra generated by those constraints, it will be very convenient to use the following 

metric per pair (/, v), 



Q 



/« _ 



E a ■ E/3 



a,/3 = l,2. 



(4.6) 



Its determinant is detQ = \E± x E 2 \ 2 , so that when Ei,E 2 are linearly independent, it has an inverse, 



1 



ffv 



|£i xE 2 \ 2 \-Ei-E- 



El 



—Ei ■ E- 
El 



(4.7) 



whose matrix elements will be denoted Q°f v - Note that this really is the analogue of the continuum 
inverse 2-metric, which we expect to appear in the Dirac algebra. It naturally induces a 'space-time 
split' as 



6> J 



l<*/8 



E> a Q%M +tA v tA v 



(4.8) 



Here we sum over a, /3 = 1,2. This identity can be used to decompose any vector A into compo- 
nents tangential to the triangle dual to v and parallel to the normal n. In the case of the flatness 
constraint/holonomy h\ v , the components are exactly the diffeomorphism and Hamiltonian constraints, 



l fv 



E i a Q'$Df ) v + n i Ht v . 



(4.9) 



The constraints are redundant in two ways: 



• Local redundancy: The holonomy has three real degrees of freedom, so that the flatness constraints 
give three constraints per face. We have define three constraints for every face-vertex pair (f,v), 
which if the two edge vectors associated to (/, v) are linear independent, are also independent. Thus 
in general the constraints associated to different vertices at the same face can be related to each 
other. 

• Global redundancy: The constraints generate vertex translations of the triangulation (which can 
be locally embedded into 3D flat space time due to the equations of motion). Given for instance a 
spherical triangulation, we can translate all vertices at once such that we either produce a global 
rotation or a global translation. Thus in this case we have six combinations among the vertex 
translations, which do not change the length and the dihedral angles of the triangulation. The latter 
provide a parametrization of the Gaufi constraint reduced phase space. Thus for the boundary of 
a tetrahedron one will have 4 (triangulation) vertices, hence 12 (flatness) constraints. Only six of 
these are independent, fixing the dihedral angles as functions of the lengths |S]. Related is the 
redundancy between the flatness constraints (on the full phase space) induced by the integrated 
Bianchi identity. For a spherical surface this identity implies that the holonomy of a given face can 
be written as a combination of the holonomies of other faces. 

The latter global redundancy has been overlooked in 29] and lead there to the conclusion that it would 
not be possible to obtain a closed (local) hypersurface deformation algebra for discrete geometries. We 
will show that this is clearly not the case. 

5 The constraint algebra 

The form of the continuum constraint algebra (I2.11[) can actually be derived by geometrical considerations 
PQ — it reflects the commutator of deformations of a hypersurface embedded in a higher dimensional 
manifold. This shows that the algebra is universal, i.e. it is independent from the field content. Moreover, 
given a certain set of assumptions on the number and metric interpretation of the (ADM) phase space 
variables, the Hamiltonian and diffeomorphism constraints can be derived uniquely |36) . (This work has 
been performed in the ADM variables, we are not aware of a derivation in connection or Ashtekar-Barbero 
variables.) 

We will also derive the constraint algebra by geometric considerations, for the boundary of a simplex 
of arbitrary dimensions, in section [TT] Let us shortly explain how the form of the algebra arises. 

(a) The commutator of diffeomorphism constraints with diffeomorphism constraints gives again diffeo- 
morphism constraints. This is due to the diffeomorphisms describing deformations in the tangent 
space to the spatial hypersurface. The commutator does not leave this hypersurface. 

(b) The commutator of a Hamiltonian and a diffeomorphism constraint gives a Hamiltonian. The 
normal to the hypersurface does not change under tangential deformations of this hypersurface. 
However a tangential vector defining the spatial diffeomorphism changes under a normal deformation 
- namely by the normal itself. Thus we obtain a Hamiltonian deformation as the commutator. 

(c) In (|5.3p we derived the variation of the normal under the change of hypersurface (in this case 
described by the fluxes). Due to the normalization of n the variation is orthogonal to the normal 
itself, that is we obtain a tangential deformation. Moreover we have the inverse spatial metric 
appearing in (I5.3[) which explains why it appears as a structure function. 

This will hold in the discrete context as well, but only for the Poisson brackets between two constraints 
for which the space-time splittings (|4.8[) are with respect to one and the same normal. This is the case for 
constraints around a single 3-valent vertex (due to the GauB constraint), as we will see in the Section [51 
However, when computing the algebra around a face, in the Section involving different vertices, hence 
different spaced-time splits, we will see that the brackets are a bit more complicated. Nevertheless, they 
remain geometrically transparent, the point being that a tangential (or normal) deformation to a triangle 
dual to a vertex v is typically not tangential (or normal) when seen in the local space-time split of a 
different triangle, say dual to the vertex v' . Therefore, a spatial deformation at v has to be non-trivially 
decomposed into diffeomorphisms and Hamiltonian at v' (and the same for normal deformations). 

Moreover, those geometric considerations only apply to the part of the algebra which is linear in the 
constraints as it refers to the geometric action that the constraints generate (we will see that the brackets 



involve terms quadratic in the constraints, which however do not generate an action on the constraint 
hypersurface) . 

We consider now the Poisson algebra between two constraints. They can in general be associated to 
two different faces /, /' and vertices v, v' . We will denote the (outgoing) edges associated to (/, v) by ei, ei 
(with indices a, (3 — 1,2) and the ones associated to (f',v') by e3,e4 (with primed indices a' , (3' = 3,4 
indices). One or both edges e3,e4 might coincide or be an inverse of one or both of the edges ei,e2. 
However, only if v = v' we will have that the fluxes associated to these edges might not commute. Also 
we need to have v C /' or v' C / (or both) for the Poisson brackets not to vanish, in addition the two 
faces need to share an edge. 

The constraints arise as contractions of the face holonomy with either the fluxes or the normal. Let 
us denote this choice by A for (/, v) and by B for (/'«')■ The Poisson brackets can then be written as 

BA l BE 1 

{A% v , B™hJ, v ,} = h l fv ^{E k , hf, v ,} B™ - h), v ,—^{E% hf v }A™ 

+ ^>h)vhy, v ,—^-^r{Ei,EP} . (5.1) 

Notice that all terms on the right hand side of this equation vanish if h\ v = h7fi v , = 0, thus the 
constraint algebra still closes (as expected). 

For A — Eg, the derivative is 8El/dE k = 5 a: p5 tk . We also need the derivative of the normal with 
respect to the edge vectors. This can be computed explicitly, but also found from the equations 

d , i dn l 

>-n) = 2n l —r = 



dE k v ' BE k a 

Cf (IT) 

■dEk^-Ep) = E ^QEk+^t = ° ■ ( 5 - 2 ) 

w a a. 

Taking the contraction of these equations with Q^E™ we conclude 

|| = -n k Q^4 . (5.3) 

The main ingredient to evaluate all the terms in (|5.1[) is the bracket between a flux at v and a 
holonomy component for (/',«'). We introduce the notation hf :vv ' for the holonomy from v to v' along 
the /' (using its orientation) and hf :v ' v for the holonomy from v' to v. Hence hf/ v i = hf :vv ihfi :v i v and 
hf v = hfi. v ' v hf. vv i. The holonomy hf v > goes along the edge e a and hence contains gj a , where Of a = ± 
is the relative orientation. The bracket between E k and hf> v > inserts the generator T k at v, the global 
sign being Of a . Thus we can write 

{E k a ,hf, v ,} = o af ,(-2)tr(h f , :vv ,T k h f ,. v , v T m )} 

= o af (-2) tr(h r ; V , v hf>: VV , T k M f ,, v , v {T m )) 

= o af , (Ad f , :v , v (T m )) n (-2) tr(h f , v T k T") . (5.4) 

Note that o a f> can be set to if e a is not in /'. Also Ad/'^'^i?) denotes the adjoint of hf :v ' v on the 
vector B (which is in the reference frame of v'), i.e. Adf>. v i v (B) = hf :v ' V B hl,. v , v . This is the transport 
of B to the reference frame of v. 

We need to be careful with the case v = v' , as then in the last line of (|5.4[) either hf>. vv ' = hf v and 
hf; V ' V = I (for o a fi = 1), or the other way around (for o a fi = — 1). In the latter case we have 

{E k ,hJ, v ,} =: o ai , (-2) tr(h f , v T k (Ad r:v , v (T m )TT n ) 

-(-2)tr(T k h f , v T m ) . (5.5) 

D— v' ,o a fr— — 1 

To evaluate the trace in the last line of (|5.4[) we expand hf v = h°f,I+h P r, v T p . We note that h a f , = 5 tr hf v 
is actually independent of the vertex v. Thus, 

{E k ,hJ, v ,} = o af ,h° f ,(Ad f ,. v , v (T m )) k + lo af ,h? flv <? kn {M f , v , v {T m )) n , 
{E k ,hf, v ,} = f o af ,h° f ,5 km +ol f ,±h p flv et> kn T n . (5.6) 



In the last line we use the notation o^t* = 1 if a C /' and o 2 a t, = otherwise. 

The evaluation of all relevant brackets, for the generic case with / ^ f and o^ti', is derived in the 
Appendix E] 

6 Constraints at a vertex 

We consider a three- valent vertex v, dual to a triangle of the triangulation. It is shared by three faces 
/)/')/"; so we have three flatness constraints. Since the triangle is closed, the flux variables meeting 
at v satisfy the GauB constraint, E\ + E 2 + E3 = 0, assuming the edges are all outgoing. A priori we 
have three different normals nf v ,nfi v ,nf" V at v. However they all coincide modulo terms proportional 
to the GauB constraint and describe the (unique) normal of the triangle dual to v. The following Poisson 
brackets are therefore modulo terms proportional to the GauB constraint. For notation and labeling of 
edges we refer to figure [5J Also we will omit the index v and instead of using an index /, /' use a prime 
to denote objects associated to the face /'. 

The simplest example of a closed triangulation with a 3-valent vertex is the boundary of a tetrahedron. 
The dual is also a tetrahedral graph, and a basis of constraints is provided by choosing three faces around 
one dual vertex v. In this way we take into account the flatness constraints of three faces. The flatness 
constraint for the last, fourth face, is redundant due to the integral version of the Bianchi identity. 
Geometrically this choice of constraint basis means that we fix one triangulation vertex in 3D space and 
translate the other three triangulation vertices normal and tangential to the triangle described by v. 
There are still global rotations left, which lead to redundancies of constraints after implementation of the 
GauB constraints. 
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Figure 2: The vertex v is 3-valent, hence surrounded by three faces, and the edges ei, e 2 , e^ meet at v. 



To find the Poisson brackets between the constraints we use the results developed in the previous 
section [3l but have to be careful to use the case v = v' of the formula (|5.6p . We find (see the Appendix 
[X]for details) for the brackets between constraints in one and the same face /, 



{D l ,D 2 } = h° f (D 1 + D 2 ) 
{H,D 1 } = -h)H + ^ det QQ 2a D a {Q w Dp 
{H, D 2 } = h° f H - ±v/det QQ la D a {Q lfi D p - 
{H,H}=0 



- Q 2p D,) 
Q 2p D?) 



(6.1) 



where the sum in a, j3 is over a, /3 = 1, 2. Also det Q = det Q a p with Q a p the inverse to Q a/3 . 

For the constraints involving two faces /, /', but one and the same three-valent vertex v = v 1 , note 
that the normals coincide modulo GauB constraints n — n' . We also have det Q — det Q' and 



(Q'Y = Q 11 + Q z - 2Q ia , (Q'r = Q -Q , (Q r r = Q 



in 



(6.2) 



The Poisson brackets between diffeomorphism constraints involving two faces /, /' are 



{Di, D' 2 } = -h°(D' 2 + D' 3 ) - \ VdelQ [Q 12 ^!!' + H(D 2 + D' 3 )) + Q 22 {D 2 H' - HD' 2 )} 
{D 1 ,D' 3 } = (6.3) 

{D 2 ,D' 2 } = (ti) a D 2 + h°D 2 . 

A Hamiltonian and a diffeomorphism constraint give 

{H, D' 2 } = h°H' - l^tQ{Q la D a ) ((Q') 2 X - (Q') 37 ^) 



{H,D' 3 } = iVdetQ (Q 2a D a ) ((Q')^D' 7 



(6.4) 



with a = 1,2 and j = 2,3. Notice that the part linear in the constraints in {H, D' 3 } vanishes. Finally 
we obtain for the brackets of two Hamiltonians from two adjacent faces 

{H,H'} = -{h') Q 2a D a -h°{Q')^D l 1 . (6.5) 

The linear part of the Poisson bracket relations mirrors the continuum algebra, i.e. two diffeomorphism 
constraints give diffeomorphism constraints, a Hamiltonian and a diffeomorphism give a diffeomorphism 
whereas two Hamiltonians give diffeomorphism constraints. Also the last Poisson bracket algebra relation 
involves structure functions in the form of the inverse two-metric. This is the general form of the Poisson 
algebra, as expected from the geometric considerations of the Section [5] and as we will see in arbitrary 
dimension in the Section [TT] This form holds due to considering constraints around a three- valent vertex, 
where the splitting into Hamiltonian and diffeomorphism constraints is with respect to the same normal 
n. When the constraints involve different normals n ^ n, additional terms in the constraint algebra will 
appear. This is the purpose of the Section [7] 

7 Constraints at one face 

We now consider the algebra of constraints on different vertices belonging to a single face (hence we 
will often drop the face label in the notations), as in the Figured] At the vertex v, we have the local 
basis (E av ,n v ) where a — 1,2 which is used to define the metric Q v and the constraints D v a = E l av h\ v , 
H" = n l h\ v . At the vertex v' , we denote the edge vectors E av i, the normal n', the metric Q v > and the 

constraints D v a , H v . 

Because we consider the brackets between constraints based on different vertices, their result is ex- 
pressed in terms of quantities which have to be transported from v to v' (or the other way around). This 
involves the holonomy hf :vv ' following the orientation of the face. We will denote E av (v') — Ad^ / w (E av ), 
n v (v') — Ad/j , (n v ). To go from v' to v, this is the holonomy hf :v / v (following the orientation of the 



face, see the Figure 3(a)), and by definition E av i(v) — Adh f . v , v (E av i), n v '(v) — Adh f . v , v (n v >). 
We find 

{D v a ,D}} = o fa h° f (Efwiv) ■ E Xv ) Q^Dl ~ o sfi hj (E av (v') ■ E Xv ,) Q$ D V J 

+ oj a h° f (nv ■ E Pv ,{v)) H v - off; h° f (n v , ■ E av {v')) H v ' . (7.1) 

The presence of H v and H v in this bracket is due to the fact that the diffeomorphim constraints D v and 
D v generate deformations in different planes (because the edge vectors (E av ) and (E av i) do not span 
the same spatial slice, as emphasized in the Figure [3(b)[ ). If the planes coincide, the prefactors in front 
of the Hamiltonians would vanish. 

The bracket between a Hamiltonian and a diffeomorphism constraint gives 

{H\D$} = -o sp h° f (rv • n v (v')) H v ' + o fa —^=e Xcr (E av ■ E 0v ,(v))DZ Q?D* 

- o f0 h° f (n v (v') ■ E Pv ,)Ql?Di - o fa h%(n v ■ Ep v ,{v))Q^ D^. (7.2) 

Among the linear part, the Hamiltonian is expected. The presence of the diffeomorphism contributions 
is again due to the fact that the local basis at v and v' do not coincide, since those terms come with the 
scalar product between the normal and the edge vectors (transported to the appropriate vertex). 
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(a) This is the face / with vertices 
V, v'. To transport a vector from 
one to the other vertex, one uses 
the holonomies following the ori- 
entation of /. 



(b) This is a gluing of triangles around a 
vertex, with the local frames associated 
to two different triangles displayed. The 
edge vectors E\jl and E^^ do not span 
the same plane, and the normals n v , n v r 
do not coincide. 



Figure 3: On the left is the face (here of degree 6) with the vertices v and v', and on the right the piece 
of triangulation (with 6 triangles) it is dual to. 

In contrast, the linear part of the bracket between two Hamiltonians only generate diffeomorphism 
constraints, 

[H\H W ] = o fa h° f (n v , ■ n v (v'))Q a J D$ - o fa h° f (n v ■ n v >(v))QfD} 

+ \o fa Q a v ?D}({n v , x n v (v')) ■ E^ V ,)Q$D V S ' - ^o fa QfD v p((n v x n v ,(v)) ■ E 1V )QZ 5 D V S . (7.3) 

Therefore the linear part of the algebra mirrors quite well the continuum algebra. In fact, the dif- 
ferences can be explained in simple geometric terms. The local bases at v and v' are different so that a 
spatial (or normal) deformation at v is typically not purely spatial (nor normal) anymore when trans- 
ported to v', but still decomposes onto D v a and H" . (Notice that the transportation plays no role, since 
one can always gauge fix, say, hf :vv / = I.) This can be understood from the view of the triangulation 
dual to the face, where it is clear that what is tangential to the triangle dual to v is not purely tangential 
to the triangle dual to v' . 

8 Geometric interpretation and constraints in Gaufi reduced 
phase space 

Holonomies and fluxes are the variables inherited from loop quantum gravity on a single graph. Since 
3D gravity deals with flat (Euclidean) geometry, we should be able to re-write the constraints and their 
algebra in terms of rotation invariant, geometric quantities, like lengths and angles, obtained from the 
holonomies and the fluxes by performing a reduction with respect to the Gaufi constraints [26] . We will 
later use these variables to generalize the constraint algebra (for the boundary of a simplex) to arbitrary 
dimensions. 

To perform the phase space reduction with respect to the Gaufi constraints we need a set of rotation 
invariant variables. They are given by the extrinsic dihedral angles (a discretized form of the extrinsic 
curvature), which are associated to the edges e, and the lengths l e of the edges e*. 




>'3 




Figure 4: Definition of the dihedral angle. 



fl 



The dihedral angle is defined as the angle between the normals of two neighbouring triangles. Here 
we need however to parallel transport one of the normals, so that we can compute the inner product in 
the same coordinate system. To avoid factor ordering ambiguities (for quantization) we choose for the 
definition of the normals the two edges which are not shared with the other triangle, see figure [4] Thus 
for the extrinsic dihedral angle at the edge e 3 in figure |4] 

maft _ (£ 4 x E 5 ) ■ g 3 {E 1 x E 2 )g^ _ N> ■ (g 3 Ng^) 

3 ~ l^x^M^x^l - \N'\\N\ (8 ' ij 



with N = Ei x Ei and N' = E4 X E5, The length l e of the dual edge e* is given as l e — y / E k Ek- 

From the Poisson brackets of the fluxes and holonomies Q3.5P one can find [26] the Poisson bracket 
relations between these variables as 

{l e ,l e ,} = , {9 e ,e e ,} = , {l e ,@e>} = See' • (8.2) 

For instance for the last Poisson bracket in the case e' — e consider the set up of figure |4j One finds that 
{E 3 -E 3l N'-(g 3 Ng^)} = 2 (g 3 (E 3 x N)g^) ■ N' ^ 2\N\\N'\\E 3 \ smQ 3 (8.3) 

from which {l 3 , <d 3 } = —1 follows. 

We will express the Hamiltonian and diffcomorphism constraints in scalar variables. This allows to 
bring the constraints into their simplest form (for the tetrahedron), which can easily be generalized to 
higher dimensions. 

Consider the Hamiltoninian H* v where v is a vertex with three outgoing edges ei,e2,e3 and / is a 
face bounded by the edges ei, e&, e 2 , see figured This Hamiltonian is given by 

(-2)tr([E 1 ,E 2 ]g^ 1 g 6 g 1 ) 



{-2)(tr(E l (g 1 g^E- 2 g & )-tT{E 2 E 1 g^g & g 1 )) . (8.4) 



\El X E; 

We will now use the flatness constraint for the face / itself in order to replace g§ with ge = g 2 g{ X ■ 
This will lead to a new constraint (H')f v , which is a combination of the constraints for the face /. It 
is however straightforward to check, that this constraint has the same action on the fluxes Ei,E 2 ,Eq 
(on the constraint hypersurface) as the constraint H^ v . In the following we will not distinguish between 
(flatness) constraints leading to same geometric action on the fluxed and drop the prime in (H')f v . This 
will not affect the terms linear in the constraints in the Poisson algebra. Thus we redefine 

H fv := lEi I ^h ((56^I56 _1 ) ■E- 2 -E 1 - E 2 ) . (8.5) 

The corresponding density two version (without the denominator) has been quantized in |19) . 
Now the exterior dihedral angle Qq can be defined as follows 

cos6 6 = N "^°^!!° } with N' = E 6 xE I , N" = E- 2 x E- 6 (8.6) 

which on the Gaufi constraint surface is equivalent to the definition (|8.1j) . Expressing the normals N', N" 
through the fluxes one can rewrite this exterior angle as 

cos p — cos a i6 cos a 2e fo _, 

cos 66 = ; : (8.7) 

smai6 sma26 

with aie, «26 the (interior) 2D angles at /* spanned between the dual edges e*, e% and eg, e 2 respectively. 
Here the angle p is defined as 

E 2 ■ {geEige 1 ) , Q „. 

cosp = — wmr ■ (8 - 8) 



3 The tetrahedron has zero physical degrees of freedom, that is the constraints determine completely the dihedral angles 
as functions of the lengths. Thus on the constraint hypersurface the action of two constraints on the dihedral angles agrees 
if the action agrees for the length variables, or equivalcntly the fluxes. 
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Compare this with the definition of the (internal) 2D angle ai2 between ef and e\ 

cosai2 = -w\m • ( } 

In the case that the 2D surface is embedded in 3D flat space the parallel transport of a vector around 
the dual vertex /* is trivial and we should have 

cosai2 = cosp . (8.10) 

Indeed this defines the 'flat' (exterior) dihedral angle as a function of the 2D angles (and hence the lengths 
of the six edges making up the three triangles meeting at /*) 

aflat, s COS a i2 - COS a i6 COS a 2 6 /qii\ 

cos6£ (ai2, ai6, a 2 6j = : ; ■ (8.11) 

smai6 sina26 

This relation generalizes to higher dimensions [37] , e.g. holds also between the 4D (exterior) dihedral 
angles and 3D (interior) dihedral angles of a 4-simplex and even for homogeneously curved simplices |38j . 
Now the Hamiltonian (I8.5[) is given as 



rjfv \ E A \ E 2 1 , *. 

H = W^W\ {cosai2 ' coap) 



(cose 6 -cos9^ at ) . (8.12) 



sma 16 sma 2 6 / f i at 

sin «i2 



The term in brackets 



C p = cos 9,3 - cos el lat (8.13) 



defines the simplest 'quadratic form' (or rather cosine form) of the constraints for a tetrahedron (with 
/3 = 1 , . . . , 6 denoting the edges) with the clear geometric meaning to impose the dihedral angles as those 
from a geometric tetrahedron. This form allows for both positive and negative dihedral angles according 
to considering the outside or the inside of the tetrahedron as the 2D surface. It appeared in [3^1 HI GO] 
and as we will see later this form can easily be generalized to higher dimensions. The 'linearized form' 
(there are two sectors ± corresponding to positive or negative orientation of the simplex) 

cf = e^ef (8.14) 

defines Abelian constraints. This follows from the Schlafli identity for the variations of the dihedral angles 
in a tetrahedron 

J2h 5e P = ° ( 8 - 15 ) 



where /? labels the edges in a tetrahedron. The Schlafli identity ensures that the dihedral angles are 
generated as derivatives with respect to the lengths from the Regge action J^b ^/3®/9 f° r one tetrahedron. 
Thus the dihedral angles are indeed the conjugated momenta to the length variables. As the dihedral 
angles can be obtained from a generating function, namely the Regge action, we have 

^ e « = w a ^ (8 - 16) 

which shows that the constraints (|8.14j) are Abelian. 

A form of the constraint (|8.12|) linear in the dihedral angles, that agrees with (|8.12j) in its action on 
the length variables on the part of the constraint hypersurface describing the positive orientation solution 
of the tetrahedron, is given by 



H = 



r\ flat • 

— smfc)g smai6Sma26 



sm ai2 
3 V 



(e 6 - @ f 6 lat ) 



- (e 6 - Q f 6 lat ) (8.17) 
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where V is the volume of the tetrahedron with edge lengths li, . . . ,l 6 and V v * is the volume of the triangle 
v* (with edge lengths l\, I2, la)- Here we fix one orientation of the tetrahedron (the boundary is given by 
the outside of the tetrahedron). Thus this constraint is not fully equivalent to the version (|8.12p . which 
allows for both orientations. Note that we can switch between the linear (in the dihedral angles) and 
cosine form of the constraints by making the exchange — sin0^' a ' -f-> cos 0. 

Similarly we can find the Hamiltonian for the face /', which is bounded by the edges e2,e4,e^~ , so 
that 



h = -iil 

h v v , 

3 V 



(©6 - S f 6 lat ) 



H' _f ' (e 4 _ e{'«*) . (8.18) 

The diffeomorphism constraints can also be deduced from their action on the length variables. One 
finds 

Di = -?i(0i-0[ iat )~; 2 cosai 2 (0 2 -0^ at )-; 6 cosai 6 (0 6 -0^ at ) 

D 2 = h{02 - Oi' at ) + h cosa 12 (0i -0{ /Qt ) + / 6 cosa 26 (0 6 -0^ a *) 

D' 2 = -/ 2 (0 2 -0^ at )-; 3 cosa 23 (03-0^ Qt )-/4COsa 2 4(04-0f at ) 

D' 3 = ; 3 (03-0^ Qt )+^ 2 cosa 23 (0 2 -0^ Qt ) + Z 4 cosa 3 4(04-0{' at ) . (8.19) 

The constraints (|8.18[) and (|8.19[) form a first class algebra. As we now have a version of the constraints 
linear in the momentum variables, this algebra will only include linear terms in the constraints. These 
reproduce the linear part of the constraint algebra found in section [6] 

9 Thiemann's rewriting 

The work [HI used a density 2 version of the Hamiltonian constraint (j8.5[) . This has the advantage that 
no factor ambiguities arise and that further more one does not need to worry of how to define the inverse 
of I Si x _E 2 |, which can be understood as the volume density of the spatial hypersurface. However to draw 
some lessons for the 4D theory one might want to address the factor ordering ambiguity as well as work 
with the density one version of the constraint, which is argued [6] to be the valid version allowing for a 
continuum limit. Furthermore these versions might differ in their action on degenerate configurations, 
which might have repercussions on the choice of boundary conditions for the physical solutions in the 
quantum theory, in particular the boundary condition for zero size hypersurfaces (the 'big bang' or 'big 
crunch'). Thus the density two version of the Hamiltonian leads to a physical wave function describing 
both orientations of the spherical universe at once [19] . 

The problem of how to divide by the inverse volume in the quantum theory was solved in 4D |6J and 
in 3D [T7] by the so-called Thiemann trick, which basically uses that 

Id/- , s 

7? = 2 d7^ M 

and expresses the derivative via Poisson brackets. These identities where derived for the continuum 
theory. We want to point out here, that such an identity can equally well be derived directly for the 
discrete geometry. 

Indeed, in the context of our discrete geometry, we can express the normalized normal of a triangle, 
spanned by dual edges e\,e\ (with ei,e 2 outgoing edges at a vertex v and the holonomy hf v starting 
with ei) aCj 

-_N l = 1e l mn R™B% where 



RT = -2tx[g^{{N-N)-^ 9l }T' 
R? = -2tr^ 1 {(N-N)i,g2}T m ) and 

N l = e\ nn E™E% . (9.2) 



4 Note that it is important to use for the normal TV* = e* m , 
via the GauB constraints. This would change the result of the Poisson bracket, which involve (71,92 and hence objects not 
invariant under rotations generated by the Gaufi constraints. 
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This gives for the contraction of Ff v (where the face / is defined by the edges e.\,e.%) with the normal 
-n* 

H^ = -Y,Fl v n k = 8tr(F ftV [R 1 ,R 2 }) (9.3) 

k 

where the square brackets are the Lie algebra brackets [R±, R2] — R1R2 — R2R1 and 

R a = Y,KT m = g-'KN-Nf^gc} (9.4) 

1 a 

with a — 1,2. This results in an expression similar to the one used in |17) for the quantum Hamiltonian. 
The difference is that here we do not necessarily average over the faces adjacent to v. For future work it 
will be interesting to compare the different factor orderings of (I9.3[) to the (density two) Hamiltonian in 

QS- 
10 4D gravity: flat sector 

The hypersurface deformation algebra of constraints can in the case of 3D discrete gravity be defined for 
arbitrary triangulated (or polygonated) 2-surfaces. It is a first class algebra - as already follows from 
the fact that the constraints arise from combinations of the flatness constraints, which are Abelian. 

However, 3D gravity is in a sense exceptional, as here discretization does not break the diffeomorphism 
symmetry of the theory - it is preserved in the form of vertex translation symmetry [321 [8] . This symmetry 
can also be translated into the quantum theory [18j [321 US] ■ This already changes if we add a cosmological 
constant, here standard Regge calculus (with flat simplices) does break diffeomorphism symmetry [15] . 
in the sense that there are no transformations that leave the action invariant and act non-trivially on 
solutions. An alternative discretization employing homogeneously curved simplices [1 51 138) can however 
be constructed, which again does preserve the symmetries. 

3D gravity is a topological theory and thus one might argue that diffeomorphism symmetry can 
only be preserved for these kind of theories. Indeed the other known examples, where diffeomorphism 
symmetry (or rather reparametrization invariance) is preserved are (0 + 1) dimensional 10 and thus 
also topological. However also in more complicated theories diffeomorphism symmetry can be preserved, 
either by admitting non-local discretizations 11 or by using the concept of cylindrical consistency to 
allow for more complicated building blocks [40] . 

Although diffeomorphism symmetry in 4D gravity is generically broken, residual gauge symmetries 
still remain [5]. For vertices with adjacent triangles carrying curvature, the vertex translation symmetry 
is broken, to an order quadratic in the deficit angle (which is proportional to the curvature). There are 
however flat solutions or even vertices which are embedded into a flat neighbourhood, for which vertex 
translation symmetry is extant. 

A canonical formulation of 4D discrete gravity can be either obtained directly from Regge calculus 
[9j [12] , or via a Gaufi constraint reduction from a loop quantum gravity like discretization, which starts 
with connection and bi-vector variables [HJ [2H1 H2] ■ 

In the analysis starting directly from Regge calculus one finds (Abelian) constraints for the linearized 
theory for every vertex, which are however changed into proper equations of motion if the non-linear 
order is taken into account [12] . Again there are special configurations, for instance a four-valent vertex 
in the 3-dimensional triangulated hypersurface, for which the constraints survive to any order. The 
reason is that such a vertex leads to a flat neighbourhood in the 4-dimensional solution. 

Indeed as was first pointed out in [39], there is a family of triangulations of the boundary of the 
4-sphere, namely boundaries of so-called stacked spheres, which lead to flat bulk solutions. For these 
triangulations we can define first class constraints. The simplest example is the boundary of the 4-simplex 
to which we will restrict in the following section. (These considerations easily generalize to constraints 
around a four-valent vertex in any triangulation.) 

The work [?TJ[35] starts, as in 3D, with a discretization using holonomies and bi-vectors. This is based 
on a discretization from the Plebanski action, which is employed in spin foams. As is well known, the 
4D case is much more challenging than 3D, due to the appearance of primary and secondary simplicity 
constraints which are mostly (in the discrete theory) second class constraints. As was first pointed out in 
[39] the reduction by these simplicity constraints can be performed in two stages. The first stage reduces 
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to a phase space analogous to loop quantum gravity restricted to the dual graph of the triangulation. 
However, this phase space is strictly bigger than the phase space corresponding to Regge calculus. The 
configurations described by this phase space were later coined twisted geometries [35]. The second stage 
involves the so-called gluing constraints 37 , which are partially second class and partially first class. 
This last first class part are known as area constraints (44j - this set is however empty for the boundary 
of a simplex. This reduction recovers the symplectic structure found directly for Regge calculus [12] 
and allows to express the constraints for the 4-simplex in the same simple form as for the tetrahedron 
(|8.13I8.14[) . Indeed the constraints are of the same form. 

Because of the complications which arise in the reduction from bi-vectors and holonomies to scalar 
variables, we will consider only scalar variables in this work, and leave the investigation of the larger phase 
space for loop quantum gravity for future work. Also, we want to point out that the Plebanski action 
agrees with the (topological) BF action to which the simplicity constraints are added. BF theory leads 
again to the flatness constraints, and the corresponding quantization on a simplex has been considered 



11 The constraints for the boundary of a simplex 

The following considerations will hold for the boundary of a rf-dimensional simplex a, with d > 3. For 
such a simplex it is convenient to introduce the following notation: We will label the vertices of the 
simplex with i = 0, . . . , d and denote by a, cr(i), cr(ij) etc. the d-simplex itself, the subsimplex of a which 
does not include the vertex i and the subsimplex of a which does not include the vertices i, j respectively. 
Correspondingly V, V(i), V(ij) will denote the volumina of these various simplices. Also 9(ij) will be 
the internal dihedral angle between the subsimplices a(i) and cr(J), hence 9(ij) is associated to the 
(d — 2)-simplex cr(ij). 

The phase space variables for the boundary of a d-simplex can be taken to be the volumin4j V(ij), j > 
i and the internal dihedral angles 9(kl), I > k: 

{V(ij),9(kl)} = 6 {iMkl) . (11.1) 

The constraints express that the dihedral angles are actually fixed as functions of the geometry: 

C {ij ) = 9 l3 9{j at (11.2) 

where 9\- ° is the geometric internal dihedral angle determined from the volumina of the simplex. These 
constraints are Abelian, again due to the Schlafli identity. Indeed in form they coincide with (|8.14j) for 3D 
gravity. More generally the same constraints and phase space variables can be defined for any ei-simplex 
(withd> 3). 

Also the geometric action of the constraints (jll.2|) is clear: it changes the geometry of the simplex 
such, that only V(ij) is altered. As we have one constraint for every subsimplex a(ij), whose volumina 
parametrize (locally) the geometry, we can change this geometry in an arbitrary way by the action of 
the constraints. Therefore with the appropriate linear combination of constraints we can reproduce the 
change of an arbitrary vertex translation. The algebra will remain first class, as we just take combinations 
of constraints. 

Thus we have to find the variation of the V(ij) under the various deformations of the geometry 
induced by vertex translations. To this end we find it convenient to introduce affine coordinates and the 
affinc metric for the simplex, which are explained in appendix |B] The dihedral angles and the various 
volumina are given by components of the inverse affine metric. 

As for the tetrahedron in the Section |H1 to obtain a set of constraints which allows for all possible 
deformations of the geometry, we can fix the vertex 0. We will consider translations of the other d vertices 
k = l,.,.,d. 

• We define d Hamiltonians H(k), one for each vertex k = 1, . . . , d, by requiring that H(k) translates 
the vertex k along the outward pointing unit normal to c(0). 



5 There are -^d(d + 1) of these volumina and -^d(d + 1) lengths variables for a cf-simplex. The volumina Viij) can be 
uniquely determined form the length variables, for d > 4 there is however a discrete ambiguity in the transition from the 
volumina to length variables. Singularities in this map appear for configurations involving orthogonal angles. Away from 
these configurations the map and its inverse can however be defined locally. 
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• The diffeomorphism constraint denoted by D(kl) will translate the vertex k along a vector that 
is tangential to ct(0) and normal to a(0l). The reason for this choice is that it corresponds to 
the geometric action of the constraints used in 4D loop quantum gravity which in form (for the 
Barbcro-Immirzi parameter that leads to a self dual connection) are the same as the 3D constraints 
(12.712.81) . 

The geometry of the flat simplex is conveniently described by its affine metric (background material 
can be found in the Appendix |Bj . The change of the affine metric fa induced by a translation of the 
vertex A: by a vector v of affine components v m is given by 

SkAfa) = diV l gij+guv l 6'? . (11.3) 

Here N(k)i = Sf is the (inward pointing) normal to the subsimplex a(i) with norm 

NWigVNWi - g kk = ^^- (11.4) 

where d is the dimension of the simplex a. It appears in (jll.31) as N(k)i(e mn f = for edge vectors e mn , 
with m,n y^ k. This ensures that the lengths of such vectors is not changed under the displacement of 
the vertex k. On the other hand we have for edge vectors e mk 

8k,v( l mk) = e lnJk,v{9ij)e 3 mk = ^ l fae° mk (11.5) 

as one expects for the change of the length square under a displacement of the vertex k by a vector v. 
Thus for the Hamiltonian H(k) we have to use a deformation vector v = —N(0) where the hat 
signifies normalization to one: 

5 H{k) {fa) = -d^-{6fS^ + ~^) =_l(^o + ^) . (n.6) 



<r 



To define the diffeomorphism constraints we need the vectors tangential to c(0) but normal to the 
subsimplex a(0l). These can be found via the induced metric for the subsimplex cr(0) 

h(p)ij = fa - fi{0)iN(p)j ■ (U-7) 

This gives the following projector onto the simplex er(0) (indices are still raised and lowered with the full 
metric) 

- ■ - ~S°g° l 

mo); = q--fr ■ ( 1L8 ) 

so that we can define the normal N(l\0), tangential to <r(0) and orthogonal to a(0l) as 

N(l\0) m = h(0) l m (11.9) 

with norm 

n lo n 10 1 WDA 2 

N(l\0) k N(l\0r = MO)" = f ~ 9 -±- = {d _ l? y\ Q y ■ (H-10) 

The diffeomorphism constraints D(kl) translate the vertex k by 

v = -N(l\G) = -(d-l)V(O) N(l\0) , (11-11) 

thus the induced change on the affine metric is given by 

fe(H)(ff«) = -(d-l)V{0)(6t% + %%-^$$ + 8ffi)\ . (11.12) 

The linear constraint Cy in (|11.2p generates changes in the area V(ij). To find the constraint 
corresponding to a deformation S k ,v we therefore have to set 

C M :=^<5 M (F(ij))C (ii) . (11.13) 



i<3 
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To find the change in the volumina V(ij), note the following identity for the volume of a simplex (see 
the Appendix IB) 

Sk.vV = ^Vg mn h,v(9mn) ■ (11-14) 

This also generalizes to subsimplices, if one uses the induced metric to project the variation of the full 
affine metric: 

6k,vV(k) = \v(k)h{k) mn 5 k , v {g mn ) 

S k ,vV(kl) = ^V(kl)h(kl) mn 5 k , v (g mn ) (11.15) 

where in analogy to (|11.7I) the induced metric on the simplex a(kl) is defined as 

h(ki) mn = h(k) mn - h{k hS )n ■ ( 1L16 ) 

Thus 

C k ,v = \ Y. V(kl)h(kl) mn tKv(grnn) (fttl ~ 0jg°*) (11.17) 

where for 8 k ,v (ffmn) we have to use (|11.6[) for the Hamiltonian constraints and (|11.12|) for the diffeomor- 
phism constraints. 

12 The simplex boundary deformation algebras 

In this way we obtain an explicit realization of the constraint algebra, however for finding the commutator 
of these constraints we do not need this representation. This is due to the linearity of the constraints in 
the momentum variables - would we use some other representation of the constraints the strategy which 
we are going to use now will only give the part of the constraint algebra that is linear in the constraints, 
i.e. the part important for the flow on the constraint hypersurface. 

To find the algebra of the constraints, we will consider the algebra of the deformations induced by 
these constraints [I], 

{/, {Ck',v',Ck,v}} — {{/, Ck>,v'},Ck,v} — {{/, Ck,v},Ck'v'} 

= [5k,v ° 5k>,v' - 5k 1 , v' ° 5k,v] (/) • (12-1) 

Knowing the action of the combination of deformations in the last line of (jl2.ll) on all metric elements 
/ = 9mn will allow us to deduce the Poisson brackets {Ck',v') Ck,v}- Applying two deformations to the 
affine metric we can write 

5k',v'(5k,v9ij) = <W [SiVj + SjvA 

= ?Wh) + ^W(vi) ■ (12.2) 

where vi — v 3 gij. Thus we have to determine how the components Vi corresponding to the Hamiltonian 
or diffcomorphism constraints change under a deformation 5. 
The easiest example is the unit normal to the simplex cr(0) 

(1 - \ 1 - fi0m~0n 

-m%) = 2 N{ ^-^ 5 ^ mn) (12 - 3) 

For the normal N(l\0)j one finds 



S(N(l\0)i) = (d-l)S[V(0) Uj 



^00 



yS) m 



m\0)i{h(Q) mn 5(g mn ) + (d-l)V(0)5]^h(0) ln 5(g mn ) . (12.4) 
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Let us consider the commutator between two Hamiltonian constraints H (k) and H(k'). For the change 
of — N(0)j under the Hamiltonian constraint H(k') we have to use (|1 1 .6|) in (|12.3|) and find 

a ok ' - 
6bw(-N{0)j) = -7=^^(0)^ ■ ( 12 - 5 ) 

This defines a deformation vector that we use in (|11.3p to get 

n ok ' /- ~ - - \ 

SH( k )oS H(kl) (~g tJ ) = |^(<^°+<^°) (12.6) 

and similarly 

-6H ik >)oS H(k) (~ gij ) = -2^(t?$ + 5?%) . (12.7) 

Comparing with 1)11.60 the sum of this terms could be interpreted.] as proportional to a Hamiltonian 
deformation at vertex k and a Hamiltonian deformation at vertex kl . The Hamiltonian and diffeomor- 
phism deformation vectors we defined are however still overcomplete. Not considering deformations at 
the vertex fixes the global translations, we have however still the global rotations left. Thus this combi- 
nation of Hamiltonian deformations can be rewritten as a combination of diffeomorphism deformations. 
To this end we have to add and subtract the appropriate term, so that 



(<W-) ° 5 H (k')(9ij) ~ S H(k') ° 5 H (k)(9ij)) = ~W " W% + Sf6f + S k 'S k 

+ #!_ (s k S° + 5 k S°) ~ ?— (S k '5° + S k 'S 
a 00 V J 3 ) a 00 V 3 3 



= (d-l)v(O) ^ r '( fefe ') ~ 6D ( k ' k )) (Sij) ■ ( 12 -8) 

Thus we obtain the Poisson bracket between two Hamiltonians 

{H(k),H(k')} = —-l——(D(k'k)-D(kk')) . (12.9) 

For the Poisson brackets between two diffeomorphisms we have to follow the same strategy and in the 
end add and subtract the appropriate term to rewrite the result into a combination of diffeomorphism 
constraints again (otherwise Hamiltonian constraints are appearing). The Poisson brackets between a 
Hamiltonian and a diffeomorphism is the simplest case, as the normal N{0)j does not change under a 
diffeomorphism. In summary we obtain the simplex boundary deformation algebra 

{H(k),H(k')} = {d _l )v{Q) (D(k'k)-D(kk')) 

{D(kl),H(k')} = (d-l)V(0)h(0) lk ' H(k) 
{D(kl),D(k'l')} = (d-l)T/(0)(/ 1 (0) fe 'Z?(fc , r)-/i(0) fc ' / '£'(fc/) + /i(0)"'(D(fcfc')-£'(A:'A : ))) . 

(12.10) 

The somewhat complicated relations involving the diffeomorphism constraints are due to the choicq3 we 
have made to use the normals to the (d — 2)-simplices <r(0l) as deformation vectors, instead of edge vec- 
tors. The appearance of structure functions for the diffeomorphisms mirrors the appearance of structure 
functions for the commutator of hypersurface normals. Thus the commutator between diffeomorphism 
constraints might provide a toy model for systems with structure functions, which could be also tested 
for the spatial diffeomorphism constraints in 3D gravity. 

On the other hand the structure function appearing for the commutator of two Hamiltonians is now 
somewhat simpler, namely just the inverse volume 1/V(0) of the 'spatial' (d — l)-simplex. 



6 Thus even the algebra involving only Hamiltonian constraints is closed. Using density two Hamiltonians should just 
add terms proportional to Hamiltonian constraints. 

7 This is motivated by the geometric interpretation of the flux variables in the 4D case. 
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Other choices for the diffeomorphism deformation vectors can be made and the algebra can be com- 
puted along the same lines. For instance we can choose, as in the Section 2] for the 3D algebra, the 
edge vectors (from vertex k to vertex k') (ekk'Y = S l k , — S k as deformation vectors. We will name the 
corresponding constraints E(kk'), which generate translations of the vertex k in the direction e^y ■ 

We find the following constraint algebra, 

{H(k),H(k')} = ^/i(0) fe '£:(fc'0-^MO) fc ''£(fcO 
i i 

{E(kl),H(k')} = (S k '-5 k k ')H(k) 

{E(kl),E(k'l')} = {Sf ~6 k ')E{kl')-(Sf,-S^)E(k'l)^dfE{kk') + 6 k E(k'k). (12.11) 

Regarding the appearance of structure functions this reflects now the continuum Poisson brackets. The 
last relation in (|12.11l) describes more generally the commutation relations between constraints generating 
translations of the vertices along the edge vectors (including 'time like' directions (eok))- Constraints 
based on such deformation vectors would lead to an algebra with structure constants instead of structure 
functions. 

So far we have been working with a set of constraints which is over-complete, even locally: at each 
vertex we have d Hamiltonian constraints from the d subsimplices of dimension (d — 1) meeting at this 
vertex. Furthermore we have dx (d — 1) diffeomorphism constraints ((d— 1) diffeomorphism constraints 
per (d— l)-dimensional subsimplex) per vertex. As a vertex can be displaced only in d directions, the set 
of Hamiltonians should be sufficient in the generic case. This will give (d+ 1) X d constraints which will 
still be a basis with global redundancy. If we subtract global rotations and global translations (or take 
into account that the length of an edge can be changed by translating both of its vertices) , we will arrive 
at \{d-\- 1) x d constraints, which is the number of edges. Thus we always have a totally constrained 
system. 

We could just choose to work with the Hamiltonian constraints only, as is done in [19]. We so far 
restricted the Hamiltonian constraints to the normal N(0), which we now relax. Thus we define the 
Hamiltonian constraints H c (kl) as deformations of the vertex k with deformation vector —c(kl)N(l). 
Here c(kl) is a normalization factor. Choosing c(kl) = 1 we obtain an Abelian algebra, as in this case the 
entries of the deformation (co-) vector —N(l)j = Sj do not depend on the geometry. Another choice for 
the normalization factors is c(kl) — d V. Thus the lengths of the deformation vector —N(l) — —c(kl)N(l) 
is the 'spatial volume' V(l). This corresponds to a density two Hamiltonian, which is the choice of [T9"] . 
Naming the corresponding constraints H(kl) the algebra is found to be 

{H(kl),H(k'l')} = dv(g kl H{k'l')-~g k ' 1 ' H{kl)\ . (12.12) 

All the algebras we have discussed in this Section hold on the boundary of a simplex in any dimension, 
starting with the two-dimensional boundary of the tetrahedroio In this sense the algebras are universal. 
We now show how they are modified in the presence of the cosmological constant. 

13 Homogeneously curved simplices 

Here we discuss gravity with a cosmological constant. We consider again the boundary of a d-simplex. 
For (2 + 1)D this can be generalized to arbitrary triangulations, while for (3 + 1)D we have again to 
restrict to the topological sector (e.g. triangulations describing stacked spheres, of which the boundary 
of a simplex is the simplest). Also we are going to discuss the theory in scalar variables, as so far a phase 
space description with (loop quantum gravity like) connection variables is not available. 

Homogeneously curved simplices provide an improved discretization of gravity with a cosmological 
constant [151 138) . In particular in 3D, whereas the discretization provided by standard Regge calculus 
breaks diffeomorphism symmetry [8j [9] , the discretization with homogeneously curved simplices is sym- 
metry preserving and triangulation independent |15} 1381 . The generalized Regge action associated to one 
homogeneously curved ci-simplex is given as (we will restrict to positive curvature, but all results can be 



8 For the boundary of a triangle we would need to associate the configuration variables to vertices. This is only possible 
if we add fields beside the metric field. 
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easily generalized to negative curvature) 






S tr = 2^V(ij)(*-0 K (ij)) + (D-l)KV . (13.1) 

(«) 

Here we describe a space of constant section curvature: the sphere with radius R — -4=. The associated 

cosmological constant is A = ' ~ '4 - - K. In (|13.1|) K (ij) is the dihedral angle at cr(ij) in the curved 
geometry as a function of the lengths of the edges. Also V and V(ij), which denote the volume of the 
simplex and of its (d — 2)-subsimplices <j{ij) respectively, have to be understood as functions of the 
lengths. The sum is over ordered pairs (mn), m < n. 

There is also a corresponding first order action [35] to (|13.ip given by 

Si" = £ WM)^ - 8(ij)) + (D- 1)kV k {6) (13.2) 

(«) 

where now the 0(ij) are independent variables and the volume V K of a is a function of the dihedral angles. 
This shows that the dihedral angles 0(ij) and the volumina V(ij) are conjugated variables. It can be 
also shown by varying the second order action (|13.1[) with respect to the length variables and using the 
Schlaefli identity 

(D-l)KSV = Y,(*j)V(ij)S9 K (tj) . (13.3) 

Hence we will adopt again the canonical Poisson brackets 

{V(ij),6(kl)} = 5 {ijUM) . (13.4) 

The constraints are similar to the flat case (see (T2] how to derive them). Like in the case k = 0, they 
relate the extrinsic curvature to the intrinsic geometry, by fixing the angles 9ij as the ones of the curved 
simplex determined by the lengths or the volumina V(ij), 

C(y) = 9ij - $ij . (13.5) 

The constraints are Abelian due to the Schlaefli identity (|13.3|) . In the curved case we can also determine 
the volumina from the dihedral angles 

C (ij) =-V(ij)+V K (ij) (13.6) 

where now V K (ij) are the volumina as functions of the dihedral angles. Those constraints are also Abelian. 
The geometric action of the constraints is that C(y) shifts V(ij), whereas C(r/) shifts the dihedral angle 
6ij. This uniquely determines the flow of the constraints on the constraint hypersurface. 

The constraints can be used to generate translations of the vertices in the homogeneously curved 
geometry, i.e. inside the sphere. We can describe those translations through deformation vectors v 
tangential to the sphere at the vertex k which chosen to be displaced. The associated constraint is given 

by 

Ck,v = y^ h,v{0{mn)) C imn) , (13.7) 

(mn) 

where the variation 5k,v{0{mn)) can be determined from the variations of the vertex vectors that we 
will introduce next. Thus for a given deformation we can find the constraints and choosing this defor- 
mation normal or tangential to the simplex boundary we obtain Hamiltonian or spatial diffeomorphism 
constraints. 

We will not need the explicit form of these constraints to compute the constraint algebra. As before 
we will instead consider the commutators of the geometric deformations induced by translations of the 
vertices. Whereas for the flat simplex we used the affine metric to encode the geometry, we now use the 
vertex vectors for a spherical simplex. The simplex is embedded in a sphere of radius -4= , which itself 

is embedded into R d+1 . Thus the vertices of this simplex can be described by (d + 1) vertex vectors ej 
pointing from the origin of M. d+1 to the vertex j — 0, . . . ,d on the sphere. We will denote by 

Gij = cos(^/nlij) = k (ej ■ ej) (13.8) 
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the length Gram matrix of the simplex. It gives access to the interior dihedral angles (and therefore the 
angle Gram matrix) as 



G ij 

\/G il Gii 



G^ = cos 9(ij) = 7=== ■ (13.9) 



Here G y is the inverse of Gj k . The matrices G,G allow to express the change of the dihedral angles 
under a variation of the vertex vectors e*, see for instance [T51 145] . 

Translating a vertex A; by a vector v just changes the corresponding vertex vector by 

$k,v(e m ) = $kmV . (13.10) 

For v to describe an allowed deformation, it has to be tangential to the sphere at e k , i.e. we have the 
condition v ■ e k = 0. 

Let us first choose deformation vectors that describe normal deformations to a given subsimplex. We 
fix the vertex j — and consider the normal to <j(0). It is orthogonal to all vertex vectors e m with m ^ 0, 
ensuring that it is orthogonal to all edges adjacent to a given vertex n =^ and not connecting to the 
vertex 0. Such a normal is also tangential to the sphere at all vertices n^O. 

These considerations suggest the dual basis (see for instance [?5]) C^W) defined by 

N(k) ■ e, = 8 kj . (13.11) 

We can then take as the (outward pointing) normal — -ZV(O) where the hat denotes normalization to one. 

We will first consider the commutator of two Hamiltonians H(k),H(k') at vertices k,k' ^ 0. This 

will suggest a choice for the vectors describing spatial diffcomorphisms. Reproducing the arguments of 

mm 

[$H(k<) °5 H (k) -^H(fe) °S H (k')] e m = 4m<5ff(fc')(-N(0)) -<5fc'm<5ff(fe)(-iV(0)), (13.12) 

we see that we have to consider the change of the deformation vector —N(0) under the translation of the 
vertices k and k' . For a general variation S acting on the dual vector N(p) we have 

S(N(p) ■ N(p)) = = 2N(p) ■ S(N(p)) 

5(N(p) -ejb) = = N(p) ■ S(e k ) + S(N(p)) ■ e k . (13.13) 



Inserting the ansatz 



S(N(p)) = aN{p) + Y J Pke k (13.14) 

k^tp 



into (|13.13[) we find that 

a = , ^ /3 m iW:P) + <5(e fc ) • 7V(:P) = . (13.15) 

m 

Here H m k(p) is the length Gram matrix for the subsimplex <r(p), i.e. excluding the vertex p. Its inverse 
has matrix elements 

H rak {p) = G mk____ (1316) 

Note the similarity to the inverse of the induced affine metric for the flat simplex (|11.10[) . Also we can 
now allow m,k = p as the inverse matrix elements are zero for these choices. 
Therefore the change of the normal N(p) is 

S(N( P )) = -Y J H km { P ) (S(e k ) ■ N(p)) e m . (13.17) 

fc,m 

Thus the normal —N(0) changes under a Hamiltonian deformation H(k) at the vertex k by 

S H(k) (~N(0)) = -J2H km (0)e m . (13.18) 
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The commutator between the two Hamiltonian deformations becomes 
[Sn(k>) $H(k) - $H(k) fe(fc')] e m = -6km ^ Hk '(°) e/ + 5k ' m X! Hkl ( Q ) e i 

= -6 km J2 Hk ' l (°) h - H ik(°) e k) + 6k> m Y, Rkl (°) ( e * - Bw (0) e v ). 
i i 

(13.19) 
In the second line, we have introduced the vectors e k i = &i — i?a-(0)efc (no summation over k). These 
vectors are in the plane spanned by e; and e k and are also orthogonal to ek- Thus the vector e k i is 
tangent to the edge between the vertices k and I at the vertex I. Consequently, it is natural to use it as 
a deformation vector at the vertex k for a diffeomorphism that we denote E(kl). Its action is 

6E(ki)e m = 8km (ei ~ H k i(0) e k ). (13.20) 

This way, we find the commutator of two Hamiltonian constraints 

{H(k),H(k')} = J2H kl (0)E(k'l)-J2H k ' l (0)E(kl) , (13.21) 

l i 

which in form agrees with the corresponding Poisson bracket (|12.11|) under the replacement of the affine 
inverse metric h kl (0) with the length Gram matrix H kl (0). Note however that as structure functions 
expressed in the volumina V(ij), they will in generally differ. This means that the algebra for curved 
simplices is a deformation of the one for flat simplices. 

Having chosen the deformation vectors corresponding to the spatial diffcomorphisms we can deter- 
mine the Poisson brackets between Hamiltonian and diffeomorphism constraints and among the diffeo- 
morphisms. We can actually generalize the vectors e k i to dki = ei — Gik&k^ where I, k can take the value 
0. The change of the deformation vector di m under a deformation 6 is 

6 (dim) = 6(e m ) - G mi 5(e l ) - k e. t (6(e m ) ■ e l + e m - 6 fa)) . (13.22) 

The resulting hypersurface deformation algebra is given by 

{H(k),H(k')} = J2H H (0)E(k'l)-J2H k ' l (0)E(kl) 
i i 

{E(kl),H(k')} = (S k '-S k 'G kl )H(k) 
{E(kl),E(k'l')} = (6f - 6 k k 'G k i)E(kl') - $, - 6 k k,Gk>v)E(k'l) - 8 k ' ' G in E(kk') + 6f,G in E(k'k) . 

Notice that in the flat limit k — > 0, Gij — > 1 and we recover the algebra for the flat simplex. However 
for generic k, we see that with our choice of deformation vectors for the spatial diffeomorphisms we 
get structure functions for the Poisson brackets between Hamiltonians and diffeomorphisms and among 
diffcomorphisms. As G k i = H k i(0) = cos( v / K/fc/) for k, I ^ the structure functions only involve the 
'spatial geometry'. 

Thus, the details of the discrete hypersurface deformation algebra does depend on the dynamics, here 
whether the cosmological constant is included or not. This differs from the continuum, where this is not 
the case. We conjecture that the algebra will in general depend on the geometry that is prescribed as a 
solution of the equations of motions. 

14 Discussion 

The Dirac's hypersurface deformation algebra can also be defined for discrete hypersurfaces. We discussed 
an explicit realization for 3D gravity in terms of connection and triad variables. Realizations of this 
algebra can also be constructed in higher dimensional gravity, if one restricts to triangulations that 
necessarily describe flat, or in the presence of a cosmological constant, homogeneously curved, geometries. 
We can draw various lessons from this work: 

• One needs to be careful to which objects to assign the Hamiltonian and diffeomorphism constraints. 
Two choices which occurred so far in the literature are (a) to associate the constraints to the vertices 
in the discretization (i.e. faces in the dual for canonical 3D gravity or 3D cells in the dual for 
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canonical 4D gravity), as is natural if the emphasis is on the action of the constraints as generators 
of vertex translations. The other option (b), used in loop quantum gravity [B] is to associate the 
constraints to the top-dimensional simplices in the spatial discretization, i.e. polygonal cells in 3D 
gravity and 3D cells in 4D gravity (corresponding to vertices in the dual) . This choice was critized 
in [35j from a discrete gravity perspective. Indeed in general the two options might constrain 
different numbers of degrees of freedom. A counting of the constraints might however be misleading 
as redundancies might occur, a point missed in [29] . 

In this work we dealt with an overcomplete basis of constraints, associated to (vertex-dual vertex) 
pairs. Starting with this choice both options (a) and (b) can be discussed using some averaging 
prescription for the constraints. A precise understanding of the underlying geometry can help to 
select the better option and to avoid cither over-constraining the theory or missing some constraints. 

• The details of the constraint algebra depend on the precise geometric definition of the constraints. 
For the Hamiltonian one has to describe the choice of a normal, for the diffeomorphism constraints 
tangential vectors to the hypersurface have to be selected. This choice influences in which relations 
structure functions will occur. 

It is possible to avoid structure functions (at least for the simplex boundary), by either using the 
Abclian deformations (|8.14|) and the one with only rescaled Hamiltonian constraints discussed above 
()12.12|1 . or deformations along the edge vectors (for the flat case). These choices do however not 
carry fully the interpretation of the continuum hypersurface deformation algebra. 

• The continuum constraints lead to a 'linear representation' of Dirac's hypersurface deformation 
algebra. We have seen that in addition, also terms of higher order in the constraints might occur. 
These terms do not lead to any flow on the constraint hypersurface. But such terms might become 
relevant if gravity is coupled to other fields. In this case one might have to reduce to a flat sector, 
even for 3D gravity. 

• The Dirac algebra is universal in the continuum, e.g. does for instance not depend on whether 
a cosmological constant is present or not. We have seen that in the presence of the cosmological 
constant the discrete constraint algebra is deformed. This is due to the discretization: the structure 
functions refer to the geometry expressed via the lengths and volumina of finite size building blocks. 
The geometry of these building blocks differs in the flat and homogeneously curved case. Thus, 
should it be possible to define constraints for discrete 4D gravity, including inhomogeneously curved 
geometries, one has to expect an even more complicated algebra. (This might be possible if the 
constraints are non-local [46J [40]). However in the limit of very small discretization scale this 
dependence should disappear, as in this case the local geometry is almost flat. 

Having clarified the classical representations for the hypersurface deformation algebra we prepared 
ground for investigating quantum representations. These seem to be in reach for the 3D theory: On 
the one hand one should investigate the algebra of quantum constraints in [19] and how this reflects 
the classical relations (|12.12|) . One should be aware that terms of higher order in the constraints could 
appear. The algebra has structure functions, thus there is the questions of how these are ordered with 
respect to the constraints on the right hand side of (|12.12|) . 

On the other hand one could consider quantization of the diffeomorphism constraints and the density 
one Hamiltonian constraints, following section [9] and adapting the techniques of [6j 117) to a simplicial 
context. Here we clarified that the Thiemann trick is also available as an exact identity in an simplicial 
phase space. Indeed, in this context one has to be careful to adjust numerical coefficients, that might 
differ from the continuum formulas. Similar considerations in 4D are also possible. 

Considering density one Hamiltonians, density two Hamiltonians or the flatness constraints will in- 
fluence how degenerate configurations are treated. This is still an open question of research [22 . This 
question is related to the choice of boundary conditions for 'zero spatial slices' and also the question 
whether one has necessarily to sum over orientations [47] . As pointed out in [19] the boundary conditions 
for the density two Hamiltonian indeed impose that the physical wave functions includes both orienta- 
tions. Thus investigating other choices for the density could uncover general mechanisms, that might also 
be applicable to bouncing cosmologies or singularity resolutions [4*8] . 

Furthermore for the quantization of the density one Hamiltonians factor ordering ambiguities arise 
that are analogous to those in the 4D theory. These might be much easier to resolve in the 3D theory, so 
that lessons can be drawn for the 4D case. 
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As we have discussed in 4D we can consider a topological sector of the discretized theory that admits 
exact constraints. These allow to define a hypersurface deformation algebra, which we here considered for 
the boundary of a simplex. Let us point out that even finding the physical wave function for the simplex 
would mean considerable progress for the 4D theory: As in the 3D theory this physical wave function 
would constitute the amplitude for the simplex form which we can built the path integral, that is spin 
foam model, of simplicial gravity. Such a physical wave function would also include information on the 
measure factors, which in 4D are not specified even for configurations describing the classically flat sector 
[49] . This of course would strengthen very much the connection between canonical loop quantum gravity 
and spin foam models, and clarify the role of diffeomorphism symmetry and the related Slavnov-Taylor 
identities along the lines of [50] . 

To this end it would be beneficial to investigate more in detail the simplicial version of the Hamiltonian 
constraints 6 , which has been started in [3T] . The question arises whether one can find a unique quantum 
representation of a given simplex boundary deformation algebra or whether there exist more than one 
representation. An answer to this question would be useful for canonical quantization as well as the path 
integral approach. 
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A Evaluation of Poisson brackets 

Here we provide some details on the generic calculation of the Poisson bracket (I5.1|) . The edges associated 
to the pair (/, v) are labeled with Greek letters a, /3 — 1,2, and those associated to (f',v') are labeled 
with a',/3' = 3,4. 

(a) Let us consider the case that A = E a . Then the first term in (I5.1[) is given by 

h l fv {El hf, v ,}B m = o af h° f ,h%(Ad f , v , v (B)) k + \o af , h p , v h% e pkn (Ad f ,.. v , v (B)) n . (A.l) 

We need to split B = (Adf/ :V / V (B)) into tangential and normal components at (f,v), 

4 



B = (B- EfiQfcE^ + (B ■ n fv )n fv . (A.2) 



Thus 

h%(Ad f , :v , v (B)) k = (B ■ EtiQftDt" + (B ■ n fv )Hf v . (A.3) 

Whereas in the first term of the right hand side of (|A.1|) we have h *, appearing which is equal 
to 1 on the constraint hypersurface, the second contribution of (|A.1|) a term quadratic in the flatness 
constraints. (Note that this term does not generate a flow on the constraint hypersurface.) To split these 
into Hamiltonian and difieomorphisms we could write 

h% = (h fv ■ E P ) Qjl E k + (h fv ■ n fv ) n% = D^Q^E* + K*« n% 
hF f , v = {h flv -E p ,)Q^El, + {h flv -n rv )n), v =D?/Qtf E%, + H?" n> f , v . (A.4) 

This would however introduce a third set of constraints at {f'v), whose splitting is not adjusted to B 
and (in the general case) not to nj v . We can transport the constraints to v' using the following rewriting 

e pkn h p f , v (Ad f ,.. v , v (B)) n = (-2) tv{h } ,.. v , v h f ,.. vv , [T k , h f ,.. v , v Bhj,] v , v }) 

= (-2) tr([B, h rv ,] (h- 1 )fv'vT k h f/v , v )) 
= e n P m B n h p f , vl (Ad {n -u, vv ,(T k )) m 

= e npm B n (Df/QfrjE*, + Hf' v 'n p f , v ,) (Ad {r) - hvv/ (T k )) m . (A.5) 
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Collecting all the terms, the first term in (|5.ip is given by (for A — E a ) 

h l fv {E l a ,hf, v ,}B m = o af , h° f , ((B • E P Q%)Df" + (B ■ n fv )Ht") 

+ \o aS , e n ? m B n (d^'q^E?, + # /V n^ v ) 
x (pFQftiMw-i^ E y ) m + H* v (Ad (/0 - 1:TO , n fv ) m ) . (A.6) 

For B equal to the normal nf v > we just need to consider the term 

e n P m n], v ,E P Y = ^—±—^(E 3 -E Y )E?-(E A -E Y )E?) 

= \E 3 xE 4 \eys'Q 5 y a v 'E™, (A.7) 

with eyg' totally antisymmetric and £34 = —£43 = 1. Using furthermore 



^,Q%Q% = det(Q^)^'= |F „; F , 2 e^' (A., 



1 
\E 3 ~E 4 \ 



(with e 34 = — e 43 = 1) wc obtain for the case A = E a and B = np v > 

h l fv {El hj, v ,}n?, v , = o af , h° f , ((n,v ■ E P Q^D^ + (n /v • n fv )H^) 

+ \o af , ] J'o'Dfr'EZ, (DJ, v Q%(Ad (n - 1:vv , E,) m + H fv (Ad (/0 -i :W n fv ) m ) . (A.9) 

Similarly we find for A = E a and B = E a i 

h l fv {E l a ,hJ, v> }E% = o af , h° f , ((E a , ■ E P QffrDf? + (E a , ■ n fv )H^) + 

§«W y/dBbQ f > v ,Qf,*e a > s , (D f /nJ, v , - H*'* ' Ef) 

x ( D i V Q^(Ad (/0 -x :w , E,) m + H fv (Ad (/0 - w n fv ) m ) . (A.10) 



(b) The other case is that A = n/„, for which we can write 

h f"jE£ i E l h ?'^ Bm = ~°«f' h °f' ®f D ^ ("/« ■ *) - 1°°/' Qfv D i V £ Pkn h P rv n) v B n . (A.ll) 

We have to rewrite the second term of the right hand side above using 

e pkn h p f , v n%B n = e^ m B n {D^Q%El+H^ , n^ vl )(M (n ^ vvl {n lv )r. (A.12) 

We find for A — Uf vi B — nf v > 

h 'fvJ^ {E k a ,hf, v ,}nf, v , = - 0af ,h° f , Qf v Dp (n fv ■ «,,„,) 



and for A — n/„, B = E a 



5°^' /, ^ ^1' D p' E y D fv qP* iA d {n -u vv ,(n fv )y n , (A.13) 



*fv-g$ {El hf, v ,}El, = -o a/ ,/»°, Q# D* (n fv ■ E a .) 



- \o aS , y/detQ f , v ,Q}Ze a , s , (d(/ nf, v , - H^'E^j D%> Q}° (Ad^-u^n/,,))™ . (A.14) 

(c) For the second term in (|5.1[) the same discussion applies. This will lead to terms linear in the 
constraints at (f'v') and to terms quadratic in constraints at (/'«') and (fv). Cancellations might occur 
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between the first two terms in (|5.1[) (for the part quadratic in the constraints if we are not considering 
/ = /' and v = v'). 

(d) The last term in (|5.1[) is only non- vanishing for v = v'. 

8vv'h l fv hy, v ,-^Qgp-{E^, E p a ,} = Y^ S wh l fv h 7'v dE k dE v i <We fcp "£" ( A -!5) 

where <5 Q7 = 1 if e a — e 7 and vanishing otherwise. This will give only terms quadratic in the constraints, 
thus cancellations might occur with the quadratic constraint terms originating from the first two terms 
in (|5.1[) . Again we can decompose the face holonomies and apply a similar strategy as in (a). 

1. A = -C/Q, B = E a , '. 

h l fv hf, v {E l a , E™,} = 6 aa ,h l fv hf, v e lmn E2 

= S aa (d^Q},^E^ + H^nJ, v ) y/detQ^Q-f feoff, (p$nj v - H*»E™) . 

(A.16) 

2. A = n fv , B = E a >: 

h l fv hf, v -g^tiE* , K,} = -Saa,, [D^Qll E% + H^nf, v ) * (Qfv)s a Q% D f *E™. 

(A.17) 

3. A = rif v , B = nf v : 

h\v h 7'v^{Et , K>}-q£t = 8aa> D$> D^ (Q fv ), a Q a / v Q a / v ~^=^ & E s ■ n f , v . (A.18) 

If we consider a three-valent vertex with outgoing edges the Gaufi constraint imposes E1+E2+E3 = 
for the three edges at v. Hence the normals to the three faces coincide and there will occur various 
simplifications, for instance case |3] above will lead to a vanishing result. 

In summary the Poisson brackets between constraints based at (fv) and (f'v') can be expressed as 
combinations of constraints at (fv) and (f'v') again. In general one has to expect also terms quadratic 
in the constraints. These terms do not lead to a flow on the constraint hypersurface, thus there is no 
geometric interpretation for these quadratic terms. The terms linear in the constraints are however 
universal and can be indeed derived from geometrical considerations as we will comment on in sections 
HI and US 

B Affine coordinates and metric for a simplex 

It is convenient to use the affine (or barycentric) coordinates introduced in [51] for Regge calculus (see 
also [52l 132] which use affine coordinates to express length derivatives of the dihedral angles and to 
proof various geometric identities) . Affine coordinates for a d-simplex are defined as follows. Let 
Vj, j = 0, . . . , d be vectors from some arbitrarily chosen point in R d to the d + 1 vertices of an d-simplex 
a. Then 

1 d 

fc=0 

define an overcomplete affine basis. Hence the affine coordinates i 3 of a vector x — Y]jX J ej are not 

— ■ 3 j 

unique. However uniqueness can be obtained by imposing the additional condition J^ ■ x 3 =0. 
A dual affine basis e 3 is defined by 

ej ■ e fc = 2} = tf - rf^i (B-2) 

where 8j is the Kronecker delta. The dual basis satisfies ^ e 3 = 0. 
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The key quantity for our calculations is the metric tensor in these coordinates. It can be shown (by 
contraction with the edge vectors) that the affine components of the metric tensor are [5T] 

k,l 
Laplace's formula for the determinant of g^ gives the squared of the d- volume of the simplex, 

V2 = J^~ ekokl ''' kd ' 1 ^ h ''' ld ' 1 ^olo---~9k d - 1 l i ^, (B.4) 

where the affine epsilon tensor is given by 

-30-u-x _ J + 1 if tnc P erm utation j Q ■ ■ ■ jd~ijd is even 
1—1 if the permutation jo ■ ■ ■ jd-ijd is odd 



and is vanishing if {jo ■ ■ ■ jd-i} can not be completed to a permutation of {0, . . . , d}. The inverse affine 



metric components g^ are defined by gu g h — g zl gik = 5 l k , and one can use (|B.4|) to write 



& - d l gifci-fcd-i xih-lt-i 5 , , g, , = L^Yl fR6) 

(d!) 3 V 2 9kd-ild-i - y 2 qj2 \ D -°) 



from which 

^=^5fe) (B.7) 

for the variation of the volume follows. 

The diagonal components of the inverse metric g u are proportional to the (n — l)-volumes of the 
(i)-subsimplices, 

9 -#^r- ( R8 ) 

The edge vectors (e mn ) 1 (from the vertex m to the vertex n) are given as (e mn ) 1 = 5 l n — 8 l m . The (inward 
pointing) normals N(i) to a subsimplex cr(i) are given as N(i)k = 5 l kl and their norm is |./V(i)| 2 = g". 
Thus one finds as a formula for the interior dihedral angles 

cos ^- TOiiiv(i)i -"7FF • ( } 
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